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We study the decay — > AtpTi'^'K~ , reconstructing the At baryon in the pK~ ■k'^ mode, using 
a data sample ol 467 x 10® BB pairs collected with the BABAR detector at the PEP-II storage 
rings at SLAC. We measure branching fractions for decays with intermediate Ec baryons to be 
1:^(2455)++^^-] =J21.3± 1.0 ± 1.0 ±5.5) x 10"^ B[B° E^{2b2Q)++pTi-] = (11.5 ± 
1.0 ± 0.5 ± 3.0) X 10~^ B[B° rc(2455)°p7r+] = (9.1 ± 0.7 ± 0.4 ± 2.4) x 10-^ and -> 
X'c(2520)'^p7r"'"] — (2.2 ± 0.7 ± 0.1 ± 0.6) x 10^^, where the uncertainties are statistical, systematic, 
and due to the uncertainty on the yl^ -> pK~Tx^ branching fraction, respectively. For decays without 
Ec{2A^^) or £'c(2520) resonances, we measure B[B° -> yl+p7r+7r-]„o„-Ec = (79±4±4±20) x 10"^ 
The total branching fraction is determined to be B[B° yl+p7r+7r"]totai = (123±5±7±32) x 10"^ 
We examine multibody mass combinations in the resonant three-particle Ecpir final states and in 
the four-particle A'^pn'^-K' final state, and observe different characteristics for the p-K combination 
in neutral versus doubly-charged Sc decays. 



PACS numbers: 13.25.Hw, 13.60.Rj, 14.20.Lq 
I. INTRODUCTION 



Decays of B mesons into final states with baryons 
account for (6.8 ± 0.6) % [1 of all B meson decays. 
Notwithstanding their significant production rate, the 
baryon production mechanism in B meson decays is 
poorly understood. Theoretical models of B meson bary- 
onic decays are currently limited to rough estimates of 
the branching fractions and basic interpretations of the 
decay mechanisms [2H6|. Additional experimental infor- 
mation may help to clarify the underlying dynamics. 

In this paper, we present a measurement of the B- 
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meson baryonic decay^ 5° — > tI+ptt+tt^. The yl+ baryon 
is observed through its decays to the pK^n^ final state. 
The study is performed using a sample of e^e~ annihi- 
lation data collected at the mass of the T(4S') resonance 
with the BaBAR detector at the SLAC National Accel- 
erator Laboratory. We include a study of the produc- 
tion of this final state through intermediate U'^^ and 
resonances. The sPlot technique [7| is used to exam- 
ine multibody mass combinations within the ScpT^ final 
states. We account for background from sources such as 
B — Dpp (mr) and B~ — > E'^pTr~ , which were not con- 
sidered in previous studies [5J 1^ . In addition, we extract 
the four-body non-resonant branching fraction and exam- 
ine two- and three-body mass combinations within the 
four-body yl+pTr+Tr^ final state. The B'^ tI+ptt+tt" 
decay has previously been studied by the CLEO [S] and 
Belle [5] Collaborations using data samples of 9.17 fb~^ 
of 357 fb~^, respectively. The present work represents 
the first study of this decay mode from BaBAR. 



The use of charge conjugate decays is implied throughout this 
paper. 
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Section [TT] provides a brief description of the BaBAR 
detector and data sample. The basic event selection pro- 
cedure is described in Sec. |III| Section [TV] presents the 
method used to extract results for channels that proceed 
via intermediate Sc baryons. The corresponding results 
for channels that do not proceed via Sc baryons are pre- 
sented in Sec. |V] Section |VI| presents the method used 
to determine signal reconstruction efficiencies, Sec. |VII| 
the branching fraction results, Sec. 
of systematic uncertainties, and Sec. 
A summary is given in Sec. IXj 



VIII the evaluation 



IX the final results. 



II. BABAR DETECTOR AND DATA SAMPLE 



The data sample used in this analysis was collected 
with the BaBAR detector at the PEP-II asymmetric- 
energy e+e" storage ring at SLAC. PEP-II operates 
with a 9GeV e~ and a 3.1 GeV e"*" beam resulting 
in a center-of-mass energy equal to the T(4S') mass 
of 10.58 GeV/ c^. The collected data sample contains 
467 X 10® BB pairs, which corresponds to an integrated 
luminosity of 426 fb~^. 

The BaBAR detector (TUj measures charged-particle 
tracks with a five-layer double-sided silicon vertex tracker 
(SVT) surrounded by a 40-layer drift chamber (DCH). 
Charged particles are identified using specific ionization 
energy measurements in the SVT and DCH, as well as 
Cherenkov radiation measurements in an internally re- 
flecting ring imaging Cherenkov detector (DIRC). These 
detectors are located within the 1.5 T magnetic field of a 
superconducting solenoid. 

Using information from the SVT, the DCH, and the 
DIRC for a particular track, the probability for a given 
particle hypothesis is calculated from likelihood ratios. 
The identification efficiency for a proton is larger than 
90% with the probability of misidentifying a kaon or pion 
as a proton between 3% and 15% depending on the mo- 
mentum. For a kaon, the identification efficiency is 90% 
with the probability of misidentifying a pion or proton 
as a kaon between 5% and 10%. The identification effi- 
ciency for a pion is larger than 95% with the probability 
of misidentifying a kaon or proton as a pion between 5% 
and 30%. 

Monte Carlo (MC) simulated events are produced with 
an e'^e~ — > BB event simulation based on the EvtGen 
program |11) and an e~^e~ — um, dd, ss, cc event simula- 
tion based on the JETSET program [12^ . Generated events 
are processed in a GEANT4 [13_ simulation of the BaBAR 
detector. MC-generated events are studied for generic 
background contributions as well as for specific signal 
and background modes. Baryonic B meson decays are 
generated assuming that their daughters are distributed 
uniformly in phase space. 



III. EVENT SELECTION 



The signal mode is reconstructed in the decay chain 
i?" A+pn+TT^ with yl+ pK~n+. All final state par- 
ticles are required to have well defined tracks in the SVT 
and DCH. Kaons and protons, as well as pions from the 
decay, are required to pass likelihood selectors based 
on information from the SVT, DCH, and DIRC. For pion 
candidates from the decay, a well reconstructed track 
is required. 

To form a A'^ candidate, the p, K~ , and 7r+ candi- 
dates are fitted to a common vertex and a probabil- 
ity greater than 0.1% is required for the vertex fit. To 
form a candidate, the A'^ candidate is constrained to 
its nominal mass value and combined with an antiproton 
and two pions with opposite charge. The mass constraint 
value differs between events from data and MC. For the 



MC events a nominal At mass of 



,MC _ 



+ 



2284.9 MeV/c2 



is chosen; this corresponds to the mass value used in 
the MC generation and to the value from fits to recon- 
structed MC events. For data, fits are performed on 
the m{pK^T:^) invariant mass distribution to find the 
nominal A^ mass. The fits are performed for each of the 
six distinct BaBAR run periods. The results are found 



to vary between m 



-.data 



data 
At 



(2285.55 ±0.18) MeV/c^ and 



(2285.62 ±0.22) MeV/c^, where the uncertain- 
ties are statistical. All invariant mass values are found to 
be consistent. The average result m'^+'^ = 2285.6 MeV/c^ 
is used as the nominal value for the mass constraint in 
data. 

Only candidates within a 25 MeV/c^ mass window cen- 
tered on the nominal Af mass m'\^}'^ 



At 



(or for sim- 

ulated events) are retained. The entire decay chain is 
refitted requiring that the direct S° daughters originate 
from a common vertex and that the probability for 
the i?° vertex fit exceeds 0.1%. 

The decays B — > Dpp (utt) with n = 1,2, which 
are described in more detail in section IIV Al can con- 



tribute a signal-like background through rearrangement 
of the final-state particles and are denoted "peaking 
background" in the following. To suppress these events, 
symmetric vetoes of ±20MeV/c^ around the nominal 
D'^ and Z)+ mass values [1] are applied in the distribu- 
tions of the invariant masses m (^[K^ 



m 



iA2 



+ \^^] B« 



and m 



where sub- 



scripts denote the mother candidate of the particles. 

To separate 5° signal events from combinatorial back- 
ground, two variables are used. The B*^ invariant 

mass is defined 



rt%. with the four- 
is ^ 



momentum vector of the B'^ candidate {Ego , P go ) mea- 
sured in the laboratory frame. The energy-substituted 
mass is defined in the laboratory frame as m-Es = 

mass energy and (£'i,Pi 



P|30 with -y/s the center-of- 
the four-momentum vector 
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of the initial e+e" system measured in the laboratory 
frame. For both variables, genuine decays are cen- 
tered at the 5° meson mass. In MC, these variables 
exhibit a negligible correlation for genuine B'^ mesons. 

To suppress combinatorial background, B'^ candidates 
are required to satisfy ttt-es S [5.272,5.285] GeV/c^. Fig- 
ure[T]shows the rriinv distribution after applying all of the 
above selection criteria. The dashed lines show sideband 



regions minv G [5.170,5.230] and 



e [5.322,5.382], 



used to study background characteristics; both sideband 
regions are combined into a single sideband region. 

The analysis is separated into two parts: I) the 



,800r 



>700- 

^600- 



|500- 
400p^,i 
300 E 
200 E 
100- 



K 



5.1 



5.15 



5.2 
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5.35 



5.4 



5.45 
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m,„„ [GeV/c^] 



FIG. 1: Distribution of the invariant mass minv [Afp-JT'^TV ) 
for events with tties in the region [5.272,5.285] GeV/c^. The 
red dotted lines indicate the signal region and the blue dashed 
lines the sideband regions. Higher multiplicity modes, such 
as B — >■ A'^pTr'^TT~Tr, appear for minv < 5.14GeV/c'^. 

measurement of the four signal decays via intermediate 
i:c(2455, 2520) resonances, i.e., B° rc(2455)++j7r-, 
B° Sc{2520)++pTr- , B° Sc{2A55fpTT+ , and B° 
I7c(2520)"p7r+, and II) the measurement of all other de- 
cays into the four-body final state A'^pTT^ir^, which are 
denoted as non-Sc signal events in the following. 



IV. B° S^°^pTr'^ ANALYSIS 

Decays via resonant intermediate states with Sc 
resonances are studied in the two-dimensional planes 
spanned by Winv and the invariant Sc candidate invari- 
ant mass m(yl+7r+) for decays with X'c(2455, 2520)++ 
and m{A+Tr-) for decays with 17^(2455,2520)°. In the 
following the like-sign AcTt invariant mass is denoted as 

and the opposite-sign invariant mass as If 

both invariant masses are referred to, we use the nota- 
tion m+±. For intermediate 172(2455, 2520)++'" states, 
B [Sc A+t:] « 100% is assumed [T]. 

We perform fits in both planes minv : w++ and 

TTiinv : 'nT'-i to extract the signal yields for the decays via 

the Sc resonances. Background contributions are vetoed 
when feasible. We distinguish between different signal 



and remaining background contributions by using sepa- 
rate probability density functions (PDF) for each signal 
and background component. We use analytical PDFs 
as well as discrete histogram PDFs. The PDFs are val- 
idated using data from the sideband regions and from 
MC samples. The different, combined PDFs are fitted to 

the minv ■ ™++ and minv ■ "^h planes and the resulting 

covariance matrices of the fits are used to calculate sPlot 
f!\ distributions of signal events. 

Figures [2ja) and [2jb) show the m_|__|_ and m^ dis- 
tributions, respectively, after applying the selection cri- 
teria as described in Sec. |III[ Signal contributions from 
the i;c(2455)++, ^^(2520)^, and Z'c(2455)° resonances 
are observed and a contribution from events with a 
^7^(2800)++ resonance is visible. The doubly-charged 
^7++ resonances are seen to contribute larger numbers 
of events than the neutral resonances. The reso- 
nant structures sit on top of combinatorial background 
and peaking background events as well as non-i7c sig- 
nal events. The latter are distributed in m(yl+7r*) like 
combinatorial background events. 



A. Background sources 

The main source for combinatorial background events 
is other B decays, while 20% originate from e+e~ — >■ cc 
events. Combinatorial events do not exhibit peaking 
structures in the distributions of the signal variables un- 
der study. In contrast, other sources of background do 
exhibit peaking structures, and are treated separately. 



1. B — > Dpp (titt) 

Decays of the type B — > Dpp (mr) with D — > (m7r), 
where n + m = 3, can have the same final state particles 
as signal decays. Rearrangement of the final state parti- 
cles can yield a fake yl+ candidate, while the 5° candi- 
date is essentially a genuine B^ suppressed only by the 
yl+ selection. Because these events represent fully recon- 
structed genuine _B-meson decays, they are distributed 
like signal events in the minv and mES variables. Table [l] 
shows the relevant decay modes and their misreconstruc- 
tion rate as signal. Furthermore, these events can also 
be misreconstructed as higher Sc resonances in the /I+tt 
invariant masses. Figure [3] shows the distributions of 
the MC-simulated background modes in the minv '■ m++ 

and minv ■ "^-| planes. Additionally, B^ — )■ D~^pp-K^ 

events with Z)+ — i4r~7r+7r+ have a minimum invariant 
mass in m++ of m (D+p) ~ 2.808 GeV/c'^ and can intro- 
duce background in the study of events with intermediate 
i7c(2800)++ resonances. 

From the misreconstruction efficiency determined from 
signal MC events and scaled with the measured branch- 
ing fractions [14], 167 ± 20 background events are ex- 
pected to contribute as signal. To suppress these events, 
veto regions are set to 20MeV/c^ around the nominal 
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FIG. 2: Event distributions in m (/I^tt^) (a) and m (/l^vr ) (b) for events in the signal region of Fig. [T] The inserts show the 
low invariant mass regions. 



D° and masses [T] in m {K~Tr^), m (ii'~7r+7r+), and 
m{K~Tr^Tr~^Tr~), with the resulting suppression rates 
given in table |Ij A systematic uncertainty is assigned to 
account for the remaining background events. No dis- 
tortions are found in other variables due to the vetoes. 
Note that D°ppn+T:^ events with D° K^Tr+ 



P 



Id" 



do 



rearranged to B — > [tt tt ' p 

fako At 

not contribute peaking background because the selection 
requirement on m {pK~ Tr~^) ^+ effectively vetoes these 
events. 



2. B° ^ (5) T^-'V+TT- 

Decays via charmonia, such as — > (g) K" 
with (5) pp and K-tt+, or (c) with 

(c) — 7> p|37r+7r~, can also produce the same final state 
particles as signal events. We observe no indication of 
such contributions in data in the relevant combinations 
of i?° daughters or in signal MC events when scaling the 
misreconstruction efficiencies with the measured branch- 
ing fractions [T]. We neglect these events, but assig n a 
corresponding systematic uncertainty (see Sec. Villi. 



as in TO++ in the ^^(2455)++ and 17^(2520)++ signal 
regions. A correlation between minv and is ap- 

parent. No significant structures are found in MC- 
generated events with nonresonant B~ — > A'^pK~Tr^ or 
with B~ Z'c(2800)+p7r~ events due to the softer mo- 
mentum constraints on the . 



4- Combinatorial background with genuine Sc events 



In both MC and data-sideband events, combinatorial 
background events with genuine i7c(2455, 2520)^+'° reso- 
nances are found to be distributed differently than purely 
combinatorial background events without E^. resonances. 
These events produce a signal-like structure in or 

m-i , but are distributed in minv similarly to purely 

combinatorial background. However, since combinato- 
rial background events with genuine i7c(2455, 2520)"'"'^'° 
resonances scale differently in minv than purely combina- 
torial background events, no simple combined PDF can 
be constructed. Thus, both combinatorial background 
sources are treated as separate background classes. 



3. B- 



S+P7T- 



5. — > A^p-K^Ti events without a Sc signal 



Events from B^ S+pir^ decays with Z'c(2455)+ 
yl+7r" or Z'c(2520)+ ^^tt" are found to have a signal- 
like shape in rriinv and ?ti_|__|_. Because of the low- 
momentum 7r° daughters in the i7c (2455, 2520)"'" center- 
of-mass systems, fake Z'c(2455, 2520)+"*" can be gener- 
ated by replacing the 7r° with a 7r+ from the S"*". Fig- 
ure |4] shows the distributions of MC-generated events. 
These events cluster in the minv signal region as well 



Events also appear as background in the m+± distri- 
bution when they contain decays into the four-body fi- 
nal state B^ — >■ A'^pTT~^TT~ , not via the signal Sc res- 
onance. For example, decays such as B^ — > S^^pn^ 
are distributed as background to B° S^pn^ events 

in m^ but as signal in m^nv Therefore, decays to 

i?° — yl+p7r'^7r~ not cascading via the signal resonance 
are included as a background class. 
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TABLE I: Efficiencies for reconstructing B" — )■ Dpp (nvr) events as signal decays by rearranging tlie final-state particles in signal- 
like combinations. In the fake signal reconstruction, the subscript particles denote the actual mother. The quantity ncxpoctod 
gives the number of fake signal events without the D-meson veto (see text), ecut gives the efficiencies of the vetoes, and ni-omaining 
gives the expected number of remaining fake events in the signal regions after applying the vetoes. The B'^ — >■ Dpp (nvr) 
branching fractions are taken from Ref. [14] and the D^/D^ branching fractions from Ref. [l]. 



Decay mode 



Fake signal 



^expected £Cut 'T'rcmaiiiing 



B° D^pp 



Bf, 



fake [P 



I BO I 



^'^^]no (6.79 ± 0.19) X 10"^ 26.0 99.3% 0.3 



D+ppTT~ 

D+ K--K+-R+ 
B° ^ D°pP!v+tv- 



5 fake ^ TT p 



Be, 



fake [TV P 



7r+] (7.28 ± 0.17) X 10-3 103.0 98.8% 1.0 



fake AZ 



fake At 



TT+l (4.19 ± 0.15) X 10-3 22.5 96.9% 0.2 



B° D*+ppK^ 
D*+ ^ D\+ 



>fake TT p 



TT+\ (2.44 ± 0.12) X 10-3 13.4 96.9% 0.1 



B. Fit Strategy 



2. Histogram PDF verification 



The signal yields of resonant decays are determined in 
binned maximum-likelihood fits to the two-dimensional 
distributions minv ■ ™++ and minv : ?7t-h Since back- 
ground events from — > Z'+pTr" decays are distributed 
similarly to signal events in all examined variables, one- 
dimensional measurements of the signal yield will not 
suffice. By extracting the signal yield in the mjnv : m^± 
plane, we exploit the fact that the distributions of ^■ 
Z7c(2455, 2520)+p7r- events are more correlated in these 
variables than signal events. 



1. Type of PDFs 

Signal and background sources are divided into 
two classes of probability density functions. Back- 
ground sources without significant correlations between 
TTiinv E^nd are described with analytical PDFs; in- 

dependent analytical PDFs are used for each of the 
two variables and a combined two-dimensional PDF is 
formed by multiplication of the one-dimensional func- 
tions. Signal and background sources with correlations 
between minv and to+± are described with binned his- 
togram PDFs Hi = Si -hi (minv, 'm+±)- For each source, a 
histogram h^ (minv,m+±) is generated from MC events, 
which takes correlations into account by design. Each 
histogram h^ is scaled with a parameter Si, which is al- 
lowed to float in the fit. Histogram PDFs are used for all 
resonant signal decays and peaking background decays 
B- 1:^(2455, 2520)+p7r-. 

In the fits to the two-dimensional distributions, the in- 
tegrals of the analytical PDFs for each bin are calculated. 
Table[lT|lists the PDFs and indicates whether they are in- 



When using a histogram PDF in fits, results prove 
to be sensitive to differences between data events and 
MC-generated events. As a cross-check, the projections 
onto TOinv are compared between data and MC simula- 
tion. The distributions are fitted using a Gaussian func- 
tion to describe signal events. The means differ between 
data and MC by A = (2.30 ± 0.25)MeV/c2. The mass 
shift does not depend on the A'^ candidate selection or 
on m^±. The most probable explanation for the differ- 
ence is an underestimation of the SVT material in the 
simulation, as studied in detail in Ref. [15j . Baryonic 
decays are especially affected by this issue, since heavier 
particles such as protons suffer more from such an un- 
derestimation compared to lighter particles. In each MC 



event, the baryon momenta 



Pp. 



and IPpgo I are there- 



eluded in the fit to mjuv : for i?° — >■ i7++p7r events 
or in the fit to m-mv ■ rn^ for — ^ S^pn'^ events. 



fore increased by 2.30MeV/c and the particle energy is 
adjusted accordingly. 

In the m4^± distributions, the means of the masses of 
the Z'c(2455)++'° baryons differ between data and MC 
by Am = (0.441 ± 0.095) MeV/c^. This effect originates 
from outdated Z'c(2455) mass inputs in the MC genera- 
tion, and so this shift is not covered by the correction for 
detector density. Z'c(2455) events are especially sensitive 
to such mass differences due to their narrow width. The 
effect is taken into account by shifting each MC event 
in m+± by +0.441 GeV/c^. The fully corrected data sets 
are used to generate the histogram PDFs employed in 
the fits to data. 



3. Combinatorial background PDF 

The combinatorial background is described by the 
PDF BGcombi Bkg term given in Table [H] It consists of 
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FIG. 3: Simulated events with 5° — > Dpp (mr) decays misidentified as signal decays in the minv : plane (left column) and 
minv :m+_ plane (right column). The MC-generated events are reconstructed as B'^ — > A^ptv'^tt' . The color scale indicates 
the relative contents of a bin compared to the maximally occupied bin (color online). 




FIG. 4: Simulated events with re(2455)+p7r" (left) and B' -> 1:0(2520) ^pvr (right) decays, where — > A'^tt^' decays 

are reconstructed as X'j""'" — >■ A'^n'^; these events accumulate in the signal regions of minv and in m++. The MC-generated 
events are reconstructed as B° — >■ A'^pTv'^-rr~ . The color scale indicates the relative contents of a bin compared to the maximally 
occupied bin (color online). 



two separable functions: for mjnv we use a first-order tion, 
Chebyshev polynomial X^Jj^^jg^^^^e ("^in^; fe) with a 
slope parameter b and for m+± a phenomenological func- 



YcombiBkg (to+±;p, g; Cup, eiow) = (1) 
(c - m+±Y ■ (m+± - eiow)' • eup . 
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The upper and lower phase-space boundaries in m+± are 
constants eiow = 2.4249 GeV/c^ and Cup = 4.215 GeV/c^. 
The phenomenological constant c = 4.108 GeV/c^ is ob- 
tained from MC and, for estimating an systematic uncer- 
tainty, varied within the values found in MC. The expo- 
nent terms p and q are allowed to float in the fits to MC 
and data. In Table ITTl S'combi Bkg is the overall scaling 
parameter of the cornbinatorial background PDF. 



4-. Combinatorial background with genuine Sc PDF 



a (m+±; Oct, b^, c^) = Ca ■ \_a„ ■ m+±^ -I- ha ■ m+± -I- l] . In 
fits to data, the width parameters are fixed to the values 
obtained from fits to MC; a systematic uncertainty on 
the shape is included by varying the parameters within 
the parameter range obtained from fits to toy MC. 

Global PDFs for fits to rriinv : and to minv : tti^ 

are formed from sums over the signal and background 
PDFs as listed in Table |TT] The global PDFs are vali- 
dated on toy MC samples of randomly chosen numbers 
of events from the signal and background classes. 



Combinatorial background events with genuine 
resonances are described by uncorrelated functions in 
TOinv and A first-order Chebyshev polynomial 

^SmbiBkgwS, (™inv;^) m mjnv is muUiplied by a non- 
relativistic Breit- Wigner function in to_|_-|- , 



Fit Results 



BWcombi Bkg w (?7l+± ; T) 
1 

(m+± - [if + 



(2) 



with mean /i, width P, and an overall scaling fac- 
tor S'combi Bkg w I to form a two-dimensional PDF 
(BGcombiBkgwE, in Table |ll| in 7ti+±. 

The PDFs for combinatorial background with and 
without genuine Yl^ resonances are validated using stud- 
ies with MC events and from fits to data within the 
rriinv sidebands of the minv : TO-i-± planes. 



fits to data distributions in 
V : planes are performed 



Maximum likelihood 
the minv : and m\ 

in the range mjnv S [5.17, 5.38] GeV/c^ and m+± S 
[2.425, 2.625] GeV/c^ covering the regions of Z'c(2455) 
and i7c(2520) resonances.^ 

The fit to rriinv : "T-H converges with ^ jvAi — 

2807/2697. Figure [5] shows the projection of the two- 
dimensional fit onto the rriinv axis. The fitted PDFs are 
shown as stacked histograms and are overlayed with the 

distribution in data. The projection onto the axis 

is shown in Fig.[6j 

In the riiinv : "r+_|_ plane, the fit converges with 
xVndf = 2592/2695. The two-dimensional fit results 
are shown in Fig. [7] for the projection onto rriinv and in 
Fig. [8] for the projection onto 



The measured signal yields are given in Table III 



Events with — >■ yl+p7r+7r^ decays but without sig- 
nal Z'c resonances are described by the product of a phe- 
nomenological function in to_|_± , 



(3) 



(rn^ 



m+±J 



and a Gaussian function Xnon-Ec (™inv; in rriinv 
Fits to mixtures of MC samples (denoted as toy MC mix- 
tures) are used to validate the combined two-dimensional 
PDF (BG^o^n"-So in Table |ll]). This procedure is designed 
to take into account the fact that — J> yl+pTr+Tr" de- 
cays without signal Ec resonances can proceed via vari- 
ous other intermediate states besides direct decays into 
the four-body final state. Since the true composition of 
B^ — )■ yl+p7r+7r~ decays without signal Sc resonances is 
unknown, toy MC sets are created by combining ran- 
domly chosen numbers of MC events with completely 
nonresonant signal decays, signal decays with intermedi- 
ate mm-Sc resonances such as N* or p, and signal decays 
with non-signal Sc resonances. 

In the fits to toy MC samples, a quadratic dependency 
on ni_|_-|- of the signal-Gaussian width in mjnv is observed. 
This is taken into account by parameterizing the width as 



D. Signal event distributions 

The distributions of signal events are extracted us- 
ing the sPlot technique [7 in variables other than 
the discrimination variables: we perform a fit to the 
two-dimensional distributions of the signal variables 
rriinv and mj^± where all shape parameters are fixed, and 
only the signal yields Ni for each signal and background 
class i are allowed to vary. Distributions for class i are 
obtained using per-event weights 



,'7i+±j = 



J2flT' Njfj (rriinv , m+± ) 



(4) 



^ Larger m+± masses are omitted because of uncertainties on 
the I7c{2800) mass values. In a study of the related decay 
B~ — > A'^p'iT~ |17| . a significant difi'erence in the mass is 
measured for the I7c(2800)° resonances with m (I7c(2800)'^) = 
(2.846 ± 0.008) GeV/c^ compared to the world averaged mass 



m (1:^(2800)0) 



PDG 



2.802 



+0.004> 
-0.007 i 



GeV/c2 [T]. Thus, the his- 



togram PDF approach based on MC simulations is not feasible 
for Sc (2800), since the input mass value is necessary for the MC 
generation. 
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TABLE II: The PDF types for signal and background sources as defined in Sec. IV B (see text for details). In the second 
column, 5*; denote scaling factors, hi histograms, and Xi, Yi, BWi analytical functions as described in the text. The third and 
fourth columns indicate in which global fit to the planes minv : rn++ or minv : rn^ a particular PDF is included. 



Mode PDF 



B" X'c(2455)++p7r" 


Hec(2455) + + — 'S's<:(2455)++ '11^^(2455)++ 


/ 




B° ^ rc(2520)++p7r- 


Heo(2520) + + ~ •S'j:^(2520)++ ' 11^^(2520)++ 


/ 




B° rc(2455)"p7r+ 


Hsc(2455)'' — 5*^^(2455)0 ' hs^(2455)0 




/ 


B° ^ X'c(2520)°p7r+ 


Hec(2520)O ~ 5*^0(2520)0 ' hj;^(2520)0 




/ 


B- rc(2455)+p7r- 


BG2o(2455)+ = 'S'^c(2455)+ ' ^£^{2455) + 


/ 




— i y (')^'?n\~^r;-Tr~ 


-L>'-iEc(2520)+ ^ '-'Sc{2520)+ 'l£c(2520)+y 


/ 

V 




Combinatorial background 


BGcombiBkg ~ ScombiSkg 

X YcombiBkg (m+±;p, q; e^p, eiow) 

, . TT-Chebvshcv / l\ 
X XcombiBkg (n^mv;6) 


/ 


/ 


Combinatorial background 


BGcombiBkgwEc = Scombi Bkg w Ec 


/ 


/ 


with genuine Sc 


X BWcombi Bkg w Ec {'m+±; fj,,r) 

.yChebyshev / 7\ 
X -^CombiBkgwEc (nT'inv'jO) 


non-^c B° -> A+pn+n- 


BGnon — Ec ~ 'S'non — Ec 

X Ynon~Ec("^+±;P, fipsb) 


/ 


/ 



X Xnj'nf Ec {rrlinv; [w^+±; Og, bg , Cg]) 




FIG. 5: Result of the fit to the minv 'Tn^ plane projected 

onto the minv axis. The data are shown as points with error 
bars; the fitted signal and background PDFs are overlaid as 
stacked histograms. 



FIG. 6: Result of the fit to the minv : rn^ plane projected 

onto the axis. The data are shown as points with error 

bars; the fitted signal and background PDFs are overlaid as 
stacked histograms. 



where Ntotai = the Ni and /i(minv ,m+±) are 

the fitted yield and PDF value for the event in the class 
i, and Vij = Cov{Ni, Nj) is the fit's covariance matrix. 
We use these weights to generate histograms in Dalitz 
variables of the ScpT^ three-body systems. Uncertainties 
are calculated for sPlot histograms with -y/^ ^ Wf . 

In Fig. |9(a)[ B° Z'c(2455)°p7r+ events are seen to 
exhibit a sharp enhancement just above the threshold in 
the m (Z'c(2455)'^7r"'") distribution; however, there is in- 
sufficient information to reliably identify the ^c(2595) + 



or ylc(2625)+ states. In Fig. 9(b) signal events from 
i?" Z'c(2455)"p7r+ accumulate in m{Ec{2Abbfp) 
only for values larger than 3.8GeV/c^, clearly ruling 
out an enhancement at baryon-antibaryon invariant 
mass threshold, which has been seen in other decays 
[m [TBHTO] . Different behavior is observed in Fig. |9(c)' 
where — Z'c(2455)°|57r+ events accumulate signifi- 
cantly only for masses m(|37r+) < 1.8GeV/c^. Possible 
sources for this structure could be A baryons or N* 
nucleon resonances decaying to p7r+. However, due to 
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mi„vfceV/c^] m(A*7i+) [OeV/c^] 



FIG. 7: Result of the fit to tlie minv -fn+j^ plane projected 
onto the minv axis. The data are shown as points with error 
bars; the fitted signal and background PDFs are overlaid as 
stacked histograms. 



FIG. 8: Result of the fit to the minv : plane projected 
onto the m++ axis. The data are shown as points with error 
bars; the fitted signal and background PDFs are overlaid as 
stacked histograms. 



TABLE III: Signal yields without efficiency correction from 
the fits to the minv : "i+± planes. The uncertainties are sta- 
tistical. 



Mode 



Signal yield 



B° ^ X'e(2455)°p7r+ 
_B° ^ X'c(2520)"p7r+ 

B° rc(2455)++p7r" 
B° X'c(2520)++p7r" 

B" -> X'c(2455)+p7r- 



347 ± 24 
87 ±27 

723 ± 32 
458 ± 38 

164 ± 104 



X'c(2520)+p7r" 273 ± 133 



are not conclusive, and are therefore not presented. 



NON-^c B 







At 



pn TV 



ANALYSIS 



or 
one- 



The rates of events decaying into the — > A^ptt^tt' 
final state without intermediate ^7^(2455, 2520)++ 
Z'c(2455, 2520)° resonances are determined in 
dimensional fits to rriinv The non-i7c signal yield compo- 
nents measured in the fits to minv '■ w++ and mjnv : rn^ 

are not used, since these yields are correlated. 



A. Fit strategy and results 



the overlap of possible broad baryon resonances, 
cannot come to a conclusion on specific modes. 

The corresponding distributions for 
I7c(2455)"'""'"p7r~ events exhibit different 



In the TO (Z'++7r ) distribution of Fig. 10(a 



behavior, 
no en- 



corresponding TO (Z'c(2455)°p) distribution of Fig. 9(b) 
making i?*^ — > Z'c(2455)"'"+|j7r~ decays more similar 



to other baryonic decays [TH [TBtil9j . In Fig. 10(c 



events from — Sc{2455)^~^p'K are distributed in 
TO(p7r~) without an obvious structure, un like events 



from -> Z'c(2455)"p7r+ in TO(p7r+) in Fig.|9(c 
The distributions from B" £'^(2550} 
in Fig. [TT] arc similar to the distributions 



from 

5° — 7> Z'c(2455)++p7r^ events shifted to higher invariant 
masses. 

Due to the relatively small event yield for 
i?" — Z'c(2520)°p7r+ decays, the corresponding sPlots 



The data are divided into two sets: 
with TO++ < 2.625 GeV/c^ and to_i 



and subset 



^11 



< 



mv 

2.625 GeV/c2, 



with 

„2 



hancement in the threshold region is visible while the 
a bump at around 2.9GeV/c^ may may be due to 
additional contributions from intermediate ylc(2880) + 
and/or ylc(2940)+ resonances. In the to (Z'c(2455)++p) 
distribution of Fig. |10(b)| events with masses below 
3.8 GeV/c^ contribute prominently, in contrast t o the 



TO++ > 2.625 GeV/c^ and 
> 2.625 GeV/c'^. Thus, potential contribu- 
tions of — > Z'c (2455, 2520) ^pTT" are confined to 



TOjjjy subset. 



In 



TO 



B° ^ Z'c(2455,2520)++'°p7r=F 



resonant signal events 
are excluded by requir- 



ing m{A+TT^) to lie outside [2.447, 2.461] GeV/c^ and 
[2.498,2.538] GeV/c^. 
In the fit of the distribution of mf^^, background from 
— > Z'c(2455)"'"p7r~ decays is taken into account using 
a double-Gaussian PDF consisting of two single-Gaussian 
functions with different means and widths. The shape 
parameters are fixed to values obtained from signal MG, 
and the signal yield is fixed to the yield measured in the 
fit to TOinv : (Table III I . Similarly, background from 
B~ — )■ Z'c(2520)+p7r~ events is described by a single- 
Gaussian PDF with fixed shape parameters from MC and 
the yield fixed to the fitted yield in TOinv : ■ Combi- 
natorial background is described with a linear function 
in TOinv Signal event contributions are described with a 
double Gaussian with a shared mean; the parameters are 
allowed to float. In the distribution of m^^^, the fit is 
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FIG. 9; Invariant mass distributions from _B° — >■ 
X'c(2455)''p7r"'" signal events extracted with the sPlot 
method. Data points are displayed in comparison with 
the distribution of reconstructed phase-space generated 
B° — >■ Z'c(2455)''p7r"'" MC events scaled to the same num- 
ber of entries (shaded histograms). 



100^ 




1.2 1.4 L6 1.8 2 2.2 2.4 2.6 2.8 3 

m(p7i") [GeV/c^] 

FIG. 10: Invariant mass distributions from 5° — >■ 
17^(2455)^^^^" signal events extracted with the sPlot 
method. Data points are displayed in comparison with 
the distribution of reconstructed phase-space generated 
B° X'e(2455)++ pn MC events scaled to the same 
number of entries (shaded histograms). 



performed with a first-order polynomial for background 
and a double Gaussian with shared mean for signal, since 
no peaking background is expected here. 



The fits are shown in Fig. 12 and the yields are given 
in Table HVI 



TABLE IV: Signal yields without efficiency correction for 
decays without intermediate X'c(2455, 2520)"*"^'" resonances. 
The first uncertainty is the statistical uncertainty. The second 
uncertainty for mf„y denotes the uncertainty on the contribu- 
tion due to -B~ — ^ EtpiT~. 



Region Signal yield 

B. Signal event distributions 

mf„^ 810 ± 88 ± 38 

Figures [13] and [14] show the combined sPlots for non- mf/^ 1918 ± 91 

Ec — )■ yl+p7r+7r~ events from the fits to m(^^ and 
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FIG. 11: Invariant mass distributions from B" — > 
£'c(2520)^^p7r~ signal events extracted with the sPlot 
method. Data points are displayed in comparison with 
the distribution of reconstructed phase-space generated 
B° -s- i:c(2520)++p7r~ MC events scaled to the same 
number of entries (shaded histograms). 



on low-mass Ec resonances appear as gaps in the distri- 
butions. We do not attempt to explicitly measure inter- 
mediate states with Z'j,(2800) baryons with the present 
approach. As described in footnote [2] significantly dif- 
fering masses of Z'c(2800) resonances have been observed 
in related B — > A^pir decays [HI [17] , which could origi- 
nate, amongst other possibilities, from different angular 
momentum states with similar masses or from contami- 
nation due to B ^ Dpp {nn) decays. Since the present 
approach uses a priori information on the masses and 
widths to generate MC-based PDFs, histogram PDFs 
cannot be applied for states with uncertain masses or 
widths. 

differ- 



In the distribution of m [pn ) in Fig. 13(c) 



ences a re seen compared to the distribution of m (p7r+) 
in Fig. 13(d) with events accumulating in m(p7r~) at 



values near the lower phase-space boundary, suggesting 
contributions from decays via A . Such a structure 
does not contribute to m{pTT~^). The r7i(7r+7r~) dis- 
tribution in Fig. 13(e) suggests an intermediate p{770) 
resonance. However, the data are not sufficiently pre- 
cise to allow a definite conclusion. The m {A'^p) dis- 
tribution in Fig. |13(f)] shows some enhancement in the 
baryon-antibaryon mass near threshold, though less 
strongly than in other measurements with baryonic final 
states, e.g., those of Ref. [T7]. A conclusive interpreta- 
tion of the m (A'^p) result is difficult, because the MC 
distribution uses all events in the allowed phase space to 
avoid a possible bias, averaging over all possible struc- 
tures. Furthermore, the projections onto the axes of the 
Dalitz space for the four-final-state-particle system make 
it difficult to identify reflections from resonances in other 
invariant masses. 

The three-body mass distributions are shown in 
Fig.[T4j Here, we do not observe structures in lower 
invariant-mass ranges that could hint at resonances, e.g., 
excited Ac baryons in m (A^ ir^ tt~ ) . The bins near the 
edges of the distribution are not reproduced correctly by 
the sPlot method, and show artificial undershoots. This 
is because the sPlot technique relies on target variables 
that are uncorrelated with the discriminating variables. 
This does not hold for points near the edges of phase 
space, where there is a dependence on TOjnv • 



VI. EFFICIENCY 



-,11 3 



In the 771 (yl+Tr"*") distribution of Fig.|l3(a) 



- . the 

contribution from intermediate Z'c(2800)++ baryons is 
clearly apparent. The corresp onding distribution in 
m {A'^'K~) is shown in Fig. 13(b) here the isospin related 
i7c(2800)° baryon is less significant. Note that the vetoes 



^ The fixed background contributions from B~ — > 
X'c(2455, 2520)+p7r~ are taken into account in the sPlot 
weight calculations, following the method described in appendix 
B of Ref. E]. 



The efficiency of the reconstruction is determined sep- 
arately for each resonant signal mode and for the non- 
Sc signal decays. Since signal MC events are generated 
uniformly in phase space, the observed decay dynamics 
are not reproduced. To avoid bias from phase-space- 
dependent reconstruction efficiencies, the MC samples 
are iteratively reweighted according to the sPlot his- 
tograms A/" [wafaJdata invariant masses for each signal 
class from the B meson daughters a, 5,. . .. The reweight- 
ing is performed iteratively over all two-daughter com- 
binations for each three-body final state and over the 
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FIG. 12: Fits to minv distributions of events decaying into the four-body final state — > ylj!"p7r+7r~ without intermediate 
X'c(2455, 2520)"'""'"''' resonances. In (a) events originate from subset mf„^ excluding the £'c(2455, 2520)"'""^'" signal regions in 
m (yl^TT*) and taking background from B~ — >■ X'c(2455, 2520)"*"p7r~ into account. In (b) the events originate from subset mf/^,. 




miAtn) [GeV/c^] m(7i;'7T;+) [GeV/c^] 




m(p7i-) [GeV/c^] m(A^p) [OeV/c^] 



FIG. 13: Two-body invariant mass distributions for B" — > A^^pn^n' signal events from the combined subsets mf^v and 
mf/^, extracted with the sPlot method. Sc resonances are vetoed in the respective invariant masses. Data points are 
displayed in comparison with the distribution of reconstructed phase-space generated — >■ Afp-jr^Tv' MC events scaled 
to the same number of entries (shaded histograms). 
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FIG. 14: Three-body invariant mass distributions for 
signal events from the combined subsets 
mfnv s-nd mf/v extracted with the sPlot method. Data 
points are displayed in comparison with the distribution 
of reconstructed phase-space generated _B° — >■ A^P'k^'k' 
MC events scaled to the same number of entries (shaded 
histograms). 



three-body combinations for the four-body non-I7c final 
states. Since we do not observe any non-trivial structures 
in the smooth and moderately varying efhciency distribu- 
tion before reweighting, we assume that the projections 
onto the overdetermined Dalitz variables are sufficient for 
reweighting. 

In the initial step i = 1, the weight is calculated 
from the sPlot histogram and the signal MC histogram 

and applied to each signal MC event. In iteration step n, 
weights are calculated accordingly from the reweighted 
signal MC from iteration n — 1 for invariant mass mbc 
with Wbc = A/'[m(,Jdata/-^["^fcc]MC^- In the following, 
the signal MC events are weighted by Wbc- Since the MC 
event values processed in step n originate from weighting 
the MC events in the previous step n — 1, the effective 
weight is w" — w"^^ ■ wtc- Negative sPlot entries are 
set to zero to avoid nonphysical weightings. 

After each step, a fit of the minv distribution of 
the reweighted signal MC is performed to obtain the 
number of reconstructed weighted events. Thus, the re- 
construction efficiency is calculated from the number of 
weighted reconstructed events and the sum of weighted 
generated events. When the reconstruction efficiencies in 
the last steps of one cycle through all B daughter com- 
binations are compatible with each other within the un- 
certainties, we assume that the reconstruction efhciency 
has converged and stop the iteration. The reconstruc- 
tion efficiencies are listed in Table |Vl The differences 
between unweighted and weighted efficiencies vary be- 
tween 2.1% for B° Z'c(2455)++p7r- and 7.9% for 
5° rc(2455)°p7r+. 



TABLE V: Reconstruction efficiencies for each signal event 
class, after reweighting signal MC events according to sPlot 
distributions. 



Mode / region 


Reconstruction cfhciency 




X'c(2455)''p7r+ 


(16.4 ±0.3)% 


5°- 


rc(2520)''p7r+ 


(16.8 ±0.3)% 




rc(2455)++p7r- 


(14.5 ±0.1)% 


B°~ 


rc(2520)++p7r- 


(17.0 ±0.2)% 






(11.6 ±0.7)% 






(16.9 ±0.1)% 



VII. BRANCHING FRACTIONS 

The product branching fractions B are calculated for 
each signal mode i with 



B 



(51 



where iVj is the sum of signal-event numbers (Tables III 
and IV I , Ei the reconstruction efficiency (Table Ivl) , and 
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the number N {T{AS)) of T{AS) decays; we assume 

that — N{r(4S)) — ~ -fo'" integrated luminos- 

ity of £ = 426fb-\ Ngg = (467.36 ±0.11 ±5.14) x 10^ 
where the uncertainties are statistical and systematic. Sc 
resonances are assumed to decay exclusively into a A^tt 
pair: B [1:^(2455, 2520)++'° ^ yl+7r±] « 100% ilj. After 
accounting for the reconstruction efficiencies, the number 
of events from the two non-Z'c(2455, 2520) decay mea- 
surements are summed and their statistical uncertainties 
are added in quadrature. 



VIII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties applying to all modes, such as 
the uncertainty on Ngg, as well as systematic uncertain- 
ties specific only to certain modes, are considered. Table 



VI 
Tot 



lists the relative systematic uncertainties Ux = 
each uncertainty x. Systematic uncertainties on the 
reconstruction efficiency of the six charged final-state 
tracks are added linearly to obtain a total tracking uncer- 
tainty. One of the largest systematic uncertainties orig- 
inates from the particle identification efficiencies. The 
uncertainties are evaluated using MC events, with correc- 
tions derived from control samples in the data. In addi- 
tion, MC events are examined without corrections. The 
relative difference in the particle identification efficien- 
cies with and without the corrections defines the uncer- 
tainty. As discussed above, Dpp (titt) decays and 
decays through charmonia states — >■ (g) K*'^ [tt+tt^] 
can yield the same combination of final-state particles 
as signal events. Based on the known branching frac- 
tions [T], a total of at most 4.5 events from these two 
event classes are expected to satisfy the signal selection 
criteria. Here, a conservative reconstruction efficiency 
of e = 0.1% is assumed, overestimating the measured 
efficiencies in signal MC. The corresponding systematic 
uncertainty {B ^ D + X, cc + X) is set equal to 100% 
of the corresponding estimated background in line 4 of 
Table |I] for each mode separately. 

Only the fits of decays via Sc resonances are af- 
fected by an uncertainty on the shape of non-Z'c B'^ — > 



A'^pn'^TT events in the mi, 



i_± planes. The param- 



eters on the signal width are varied individually by one 
standard deviation and the maximum deviation in the 
event yield is taken as the systematic uncertainty (la- 
beled "nonres. shape" in Table VI ) . For the shape of com- 



binatorial background, a systematic uncertainty is de- 
termined by varying in the PDF the constant describ- 
ing the end point of the phase space, Cup- Fits are re- 
peated with the constant moved from the nominal up- 
per phase-space limit towards the upper end of the fit 
region in 0.2GeV/c^ steps. The maximal deviations in 
the fitted signal event yields are taken as the systematic 
uncertainty, labeled "Combi Bkg shape" . For histogram 
PDFs, systematic uncertainties are negligible due to the 
large size of the input MC data sets. In studies on con- 



trol variables, we found a good agreement between data 
and MC-generated events. For deviations of the MC- 
generated events from data, we applied corrections as 
described in |IVB 2| An uncertainty, labeled "Eff. Corr.", 
on the efficiency-calculation weighting is evaluated after 
completing a cycle through all daughter combinations for 
each mode. The values converge and are within the sta- 
tistical uncertainties for all modes after one full cycle, 
except for 5° Z'c(2520)°p7r+. For this mode, the ef- 
ficiencies differ by 1.9% with a statistical uncertainty of 
1.8% and the difference is taken as additional systematic 
uncertainty. 

A systematic uncertainty, labeled "i?^ (/ + //)", on 
the contribution of B^ — > Z'c(2455, 2520)+p7r^ events 
to the non-Z'c — yl+pTr+Tr" yields is calculated by 
repeating fits to the m(^^ distribution assuming no con- 
tribution and overestimating the found contribution by 
a factor of two. The maximum deviation of 38 events is 
included as systematic uncertainty. 



IX. RESULTS 



Table |VII| lists the measured branching fractions 
for each mode. Due to the large uncertainty on 
B [A'^ ^ pK^?:^], its systematic uncertainty is given 
separately. The uncertainties on the total branching frac- 
tion of all B° — >■ ^^pTT+TT^jg^j decays are added quadrat- 
ically when uncorrelated and linearly when correlated. 

Because the — )■ I7c(2520)'^p7r+ signal has less than 
three standard deviations significance, we also report 
a 90% confidence level (CL) upper limit for this chan- 
nel. The upper limit is determined using Bayesian meth- 
ods, with statistical and systematic uncertainties added 
in quadrature. We do not include the uncertainty on 
B[A^ pK^TT^] in the systematic uncertainty of the 
upper limit, but factor out the current branching frac- 
tion of 0.05 H]. Assuming a Gaussian distribution, the 
90% integral of the physically meaningful region B > 
yields 



B[B'' i;e(2520) 



B[A+ ^ pK-n+] 
005 
< 3.10 X 10"^ 



(6) 



at 90% confidence level. 

Resonant decays via I7c(2455, 2520) baryons provide a 
large contribution to the four-body final state. The ratios 
of decays via Sc resonances in comparison to the largest 
such mode, B° Z'c(2455)++p7r~, are 



B [B° rc(2455)V+] 
B [BO -> i:^(2455)++p7r- 

B [B° Z'c(2520)++p7r- 
B [B° Z'c(2455)++p7r- 



= 0.425 ± 0.036, (7) 



0.541 ±0.052, (8) 
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TABLE VI: Summary of relative systematic uncertainties Ui = in [%] for non-X'c — > A'^pn^n , resonant — >■ 

X'c(2455)p7r*, and resonant — >■ X'c(2520)p7r* decays. The total systematic uncertainties, added in quadrature, are given in 
the last row. 



Uncertctinties Ui [9o] 




£'c(2455)" 


X'c(2455)++ 


X'c(2520)" 


£■£(2520)++ 




1.1 


1.1 


1.1 


1.1 


1.1 


Tracking 


1.2 


1.2 


1.2 


1.2 


1.2 


PID 


4.3 


4.3 


4.3 


4.3 


4.3 


B ^ D + X, CC + X 


0.1 


0.5 


0.5 


0.5 


0.5 


nonres. shape 




0.1 


0.1 


0.4 


0.4 


Combi Bkg shape 




0.001 


0.001 


0.7 


0.7 


Eff. Corr. 








0.1 




B- (7 + //) 


1.8 












4.9 


4.6 


4.6 


4.7 


4.7 



while the fraction of all decays that proceed via the 
Z'c(2455)+p7r^ mode is 



B [B° Z'c(2455)++p7r-] 

total 



= 0.174 ±0.047. (9) 



In these three results, systematic uncertainties common 
to numerator and denominator cancel, and only the sys- 
tematic uncertainties specific to each mode are added 
in quadrature. The three-body intermediate states have 
comparable branching fractions to the non-resonant three 
body decays B A+pir |H1 [TCI [HI HO] ■ 

The measured branching fractions are in good agree- 
ment vi^ith previous measurements from Belle [9j. 



X. SUMMARY 



We observe the decay B° — )• A^pir^ir , study 
the intermediate decays via 17^(2455)++, 17^(2520)++, 
X'c(2455)'^, and Z'c(2520)° resonances, and measure 
their branching fractions (Table [VlT ). Yields for events 
decaying to — ?> tI+ptt+tt" without intermediate 
l^c(2455, 2520)'*'''"^° resonances are obtained from one- 
dimensional fits in minv, taking information from fits 
to minv : fn+± into account. For all decay modes, we 
show the sPlot distributions of the signal, and we observe 
significant differences between decays into S'^^p'K~ and 



S^piT^ final states. 
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